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Quantum systems are ultimately characterized by the ob- 
servable correlations they exhibit. For instance, an observ- 
able such as the correlation function between two spins in a 
typical spin chain may decay exponentially as a function of 
the distance separating them or, in the case the system under- 
goes a phase transition, algebraically. The coiTect assessment 
of these quantum correlations is tantamount to understanding 
how entanglement is distributed in the state of the system. 
This is easily understood as follows. Let us consider a con- 
nected correlation 



(vl/|0,0,|*),^ 

(*|o.o,|*) 



(1) 



(vl'|0,|vl/)(*|0,|vl') 



where Oi and Oj are operators at sites i and j respectively. 
This connected correlator would vanish identically for any 
product state l^*) = That is, Oi ®Oj is a product op- 

erator and, consequently, its coiTelations can only come from 
the amount of entanglement in the state \^). It follows that 
the ground state of any interesting system will be highly cor- 
related and, as a particular case, even the vacuum displays a 
non-trivial entanglement structure in quantum field theories. 

Notice that, at this point, our emphasis has moved from 
Hamiltonians to states. It is perfectly sensible to analyse en- 
tanglement properties of specific states per se, which may be 
artificially created using a post-selection mechanism or may 
effectively be obtained in different ways using various inter- 
actions. We are, thus, concerned with the entanglement prop- 
erties that characterize a quantum state. Yet, we shall focus on 
states that are physically relevant. In particular, we shall study 
the entanglement properties of ground states of Hamiltonians 
that describe the interaction present in spin chains. 

It is clear that the property of entanglement can be made 
apparent by studying correlations functions on a given state. 
We could consider two-, three- or n-point connected correla- 
tion functions. Any of them would manifest how the original 
interactions in the Hamiltonian have operated in the system to 
achieve the observed degree of entanglement. For instance, 
free particles (Gaussian Hamiltonians) produce n-point cor- 
relators that reduce to products of two-point correlators via 
Wick's theorem. Nonetheless, the study of specific correla- 
tion functions is model dependent. How can we compare the 
correlations of a Heisenberg Model with those of Quantum 
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Chromodynamics? Each theory brings its own set of local 
and non-local operators that close an Operator Product Expan- 
sion. Different theories will carry different sets of operators, 
so that a naive comparison is hopeless. A wonderful possibil- 
ity to quantify degrees of entanglement for unrelated theories 
emerges from the use of Renormalization Group ideas and the 
study of universal properties. For instance, a system may dis- 
play exponential decays in its correlations functions which is 
globally controlled by a common correlation length. A model 
with a larger correlation length is expected to present stronger 
long-distance quantum correlations. 

We may as well try to find a universal unique figure of merit 
that would allow for a fair comparison of the entanglement 
present in e. g. the ground state of all possible theories. Such 
a figure of merit cannot be attached to the correlations proper- 
ties of model-dependent operators since it would not allow for 
comparison among different theories. The way to overcome 
this problem is to look for an operator which is defined in ev- 
ery theory. It turns out that there is only one such operator; the 
stress tensor. To be more precise, we can use the language of 
conformal field theory which establishes that there always is 
a highest weight operator that we call the Identity. The Iden- 
tity will bring a tower of descendants, the stress tensor being 
its first representative. Indeed, the stress tensor is always de- 
fined in any theory since it corresponds to the operator that 
measures the energy content of the system and it is the oper- 
ator that couples the system to gravity. Correlators of stress 
tensor operators are naturally related to entanglement. In par- 
ticular, the coefficient of the two-point stress tensor correlator 
in a conformal field theory in two dimensions corresponds to 
the central charge of the theory. 

There is a second option to measure entanglement in a given 
state with a single measure of entanglement which is closer in 
spirit to the ideas of Quantum Information. The basic idea 
consists of using the von Neumann (entanglement) entropy of 
the reduced density matrix of a sub-part of the system which 
is analysed. Indeed, the entanglement entropy quantifies the 
amount of surprise that a sub-part of a system finds when dis- 
covering it is correlated to the rest of the system. Therefore, 
entanglement entropy is a bona fide measure of the correla- 
tions in the system. The advantage of the von Neumann en- 
tropy of entanglement is that it can be defined for any system. 
We expect its general properties, as the way it scales with the 
size of the sub-part of the system we are considering, should 
characterized the quantum state in a quite refined way. 

It is tantalizing to exhaustively explore the behaviour of the 
entropy of entanglement in relevant physical systems. For in- 
stance, will the entropy of entanglement scale differently at a 
critical point as compare to a non-critical one? Will scaling 
properties depend on the dimensionality of the system. Is dis- 
order relevant for long-distance correlations? Are there non- 
local systems where entropy obeys some singular behaviour? 
How does entanglement renormalize? How does entangle- 
ment evolve dynamically? We can even go further away from 
standard dynamical models and question whether entangle- 
ment is somehow related to computational complexity prob- 
lems, both NP-complete and QMA-complete. We shall now 
briefly review some of these questions. 



II. AN EXPLICIT COMPUTATION OF ENTANGLEMENT 
ENTROPY 

Let us start our discussion with the study of the behaviour 
of entanglement at different regimes (critical and non-critical) 
of the XX model. As we shall see, entanglement entropy will 
be a good tool to describe the properties of the quantum phase 
transition which characterize this model 

In order to do this, first we need to introduce the Von Neu- 
mann entropy as a measure of the bipartite entanglement in 
pure sates. Then, we will study the scaling of the entangle- 
ment entropy for the simple XX model. We then will proceed 
to compute the ground state \GS) of the system, from which 
we can obtain the spectra of the reduced density matrix pL for 
the block of L contiguous spins. The knowledge of the eigen- 
values of PL will let us determine its entanglement entropy 
Sl- Finally, we are going to analyse how the entanglement 
behaves depending on the critical properties of the model. 



A. Entanglement Entropy 

The problem of measuring and quantifying quantum corre- 
lations, or entanglement, in many body quantum systems is a 
field of research in its own, that benefits both from condensed 
matter and quantum information ideas. Here, we shall only 
discuss the Von Neumann entropy as a figure of merit for en- 
tanglement. Nevertheless, there are many other measures that 
have been largely explored. A detailed explanation of them 
can be found in several reviews 

Our choice for the entropy of entanglement is based on a 
combination of ideas. Entropy has a clear information the- 
ory meaning. It also relates to extensive research in quantum 
field theory and the physics of black holes. Furthermore, its 
scaling properties are related to the characterization of quan- 
tum phase transitions as provided by conformal field theory. 
On the other hand, entanglement entropy is not a simple quan- 
tity to compute, neither a direct observable (though it relates 
to them). 

The von Neumann entropy can be used to measure the en- 
tanglement between two parts of the system, that we call A 
and B. Let us take a ket ab belonging to 7i = ® Hb- 
According to the Schmidt decomposition, for any pure bipar- 
tite state we can always find two orthonormal basis {|<;5i)yi} 
and {|(/)j)b} such that the state \iP)ab can be written as 

X 

\tpAB) ='^ai\(p^)A\(pz)B , (2) 
i 

where can be chosen real and positive and are called 
Schmidt coefficients, and x < min (dimTi^ , dimTi s ) is the 
Schmidt number. Note that the Schmidt decomposition is just 
the diagonalization of the matrix of coefficients in the original 
state which is always possible if we can perform two indepen- 
dent unitary transformations in A and B. 

The Von Neumann entropy between these bipartitions is de- 
fined as the Shannon entropy of the square of the Schmidt co- 
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efficients. 



Sa = Sb ^ - ^^Q^ log 



(3) 



This expression can be written in terms of the reduced den- 
sity matrices of each part of the system. That is, 



Sa = S{pa) = -tr (pa log2 Pa) , 



where 



PA = tr Bi\ip) 



AB\ 



\AB 



'a1\(j)i)B{'>Pi 



B ■ 



(4) 



(5) 



It is easy to see that Sa = Sb- Thus, the surprise that A 
experiences when discovering its correlation to B is identical 
to the one of B realizing its correlation with A. 

The von Neumann entropy verifies the following properties: 
/) it is invariant under local unitary operations {Sa ~ Sb is 
a function of the a^'s only); ii) it is continuous (in a certain 
sense also in the asymptotic limit of infinite copies of the state; 
see e.g. Ref. 01); //;) it is additive: S{\i}j) ® |0)) = S{\^j)) + 

sm- 

In our particular case, we are going to use the entanglement 
entropy to study the quantum correlations of spin chains. We 
will be interested in determine the entanglement between a 
block of L contiguous spins and the rest of the chain. Then, 
if \GS) represents the ground state of a system of N spins, 
PL = tr n^l{\GS) {GS\) is the reduced density matrix of 
the block of L contiguous spins that we will use in Eq. |4] 

Finally, let us point out that, in the case the ground state 
that we study is translationally invariant, neither p^ nor Sl 
will depend on the position of the block of spins in the chain. 
In this case, it is easy to show that the entropy Sl is a concave 
function respect to L lfl2ll 



Sl > 



Sl-m + Sl+m 



(6) 



where L = 0, • • • , iV, and Af = 0, • • • , min {N-L, L}. 



B. XX model 

We shall now present a computation of entanglement en- 
tropy for the reduced density matrix of the ground state of the 
widely studied XX model llll^l- This theory captures the non- 
trivial structure of a quantum phase transition, while remain- 
ing simple enough to carry explicit computations throughout. 
The XX model consists of a chain of N spin i particles with 
nearest neighbour interactions and an external magnetic field. 
Its Hamiltonian is given by 



N~l 



N-1 



Hxx = - J E i'^f'^f+i + ^^f+i) + E ' (7) 



where I labels the N spins, A is the magnetic field and 
(/i = X, y, z) are the Pauli matrices at site I. 



Without loss of generality, we are going to consider that the 
magnetic field is oriented in the positive Z direction (A > 0), 
since, if this was not the case, we could always map the system 
onto an equivalent one with A > by simply interchanging 
the spin states up and down. 



C. Ground State 

Next, we need to compute the ground state \GS) of the XX 
Hamiltonian (|7|. In order to do this, we will follow two steps: 
(/) first, we will perform a Jordan- Wigner transformation to 
rewrite Hxx as a quadratic form of fermionic operators, and 
then (//) we will take profit of the translational invariance of 
the system realizing a Fourier transform which will diagonal- 
ize the Hamiltonian. A third step which is needed in the more 
general XY model, the Bogoliubov transformation, is not nec- 
essary in this particular case. Let us remark that this cornpu- 
tation is standard and appears in many text books iflslfTilflsll . 

The Jordan- Wigner transformation maps a spin chain of 
interacting qubits onto an equivalent system of interacting 
fermions. This powerful transform is defined by the follow- 
ing relation between the Pauli matrices and the creation and 
annihilation of the fermionic modes 



ai 



' i-i 

n 



(8) 



We, indeed, can check that the fermionic operators ai fulfil 
the canonical commutation relations 



{a\, am} = Shn, {ai, a,n} = . 



(9) 



The idea behind the transformation is to identify the state of 
the spin Z (0 or 1 in the computational basis) with the occu- 
pation number of the corresponding fermionic mode. Thus, 
in Eq.[8j the factor {af — iaf)/2 corresponds to the operator 
|0) (1| in the computational basis, and the product Hm^o '^m 
generates the appropriate sign in order to satisfy the commu- 
tation relations. 

The Jordan- Wigner transformation casts the XX Hamilto- 
nian onto 



Af-l 



Af-1 



Hxx = - E + o-l+io-i) + E '^i"' ' ^^^^ 



l=Q 



1=0 



which corresponds to a model of free fermions with chemical 
potential A. 

Now, let us exploit the translational symmetry of the system 
by introducing the Fourier transformed fermionic operators 



1 JV-l 

6. . E 



(11) 



where < k < N — 1. As the Fourier transform is a unitary 
transformation, these new b^. operators also satisfy the canon- 
ical commutation relations and, therefore, they are fermionic 
operators. 
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The Hamiltonian, written in terms of these operators, 
displays a diagonal structure 



Hxx = J2 ^'^^i^" 



(12) 



fc=0 



where the energy that penalizes (or favours, depending on the 
sign) the occupation of mode k is 



At = A — 2 cos 



27rfc 



(13) 



We have assumed that the system satisfied periodic boundary 
conditions. If this was not the case, the Hamiltonian would 
not be diagonal due to an extra term proportional to -i-. In the 
thermodynamic limit, therefore, this extra term disappears. 

We realize that, on one hand, if A > 2, then > V fc. 
This implies that the ground state of the system is the state 
annihilated by all bk operators 



bk\GS) = Vfc, 



(14) 



and, therefore, it has energy. 

On the other hand, if 2 > A > 0, the ground state is the 
state annihilated by the operators bk with Afc > and bl„ with 

A,„ < 0, 



bk\GS)=0 if Afc > 
bl\GS)=0 if A,„ < 0, 



(15) 



"AT 




— arccos 


m 


_2tt 





and its energy is simply J2m ^™ ^ ^™ < ^- Fig.[T]and Eq. 
( fT3] l. we can see that if fee > ^ > or A^ - 1 > fc > A^ - fee, 
where kc is defined by 



(16) 



then Afc < 0, whereas for the rest of cases Afc > 0. In Eq. 
( fT6] l, the brackets [] represent the floor function. 
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FIG. 1: The two terms of A^, Eq. il3\ . are plotted for the particular 
case A = 1. We realize that if 2 cos (2^) > A, Afc < while if 

' 27rfe ~l 



2 COS (2^) < A,Afc >0. 



D. Entanglement entropy of a block 

The strategy to get the Von Neumann entropy of a block 
of L spins first consists in computing the correlation matrix 
(aj„a„) of the GS in this block. Then, the eigenvalues of this 
correlation matrix are related with the eigenvalues of the re- 
duced density matrix of the block which are required to deter- 
mine the entanglement entropy. 

The simple structure of the GS, shown in Eqs. (fl4] i and ( fTSl ). 
makes easy to compute its correlation matrix 



if Ap < 
if Ap > 



(17) 



From now on, we will consider the case in which 2 > A > 0. 
Notice that if A > 2, then (6^65) = for all p and q. This 
case is trivial to analyse, since the correlators (aj„a„) are also 
null, and the GS is in a product state. 

The next step is to go back in the Fourier transform to get 
the correlation matrix of the a„ operators 



{aLan) ]^ 2^ cos 

fe=0 



27r, 



fc(m — n) 



(18) 



In the thermodynamic limit, the previous sum becomes an in- 
tegral and it can be determined analytically. In this case, the 
correlation matrix of the block of L fermions in position space 
is, 



{oLan) 



1 sinfcc(TO 

TT TO — 7 



■n) 



(19) 



where L > m,n > 0. Notice that, by means of Wick's the- 
orem, any operator that acts on the block can be written in 
terms of the correlation matrix Amn- For instance, 



(20) 



This is due to the fact that the system is Gaussian, and its 
eigenstates are determined by the first and second moments of 
some fermionic operators. 

The correlation matrix Amn could also be computed using 
the density matrix of the block p^. 



Ari 



Tr{ajnanPL) ■ 



(21) 



We, thus, need to invert the previous equation, that is, to com- 
pute the density matrix from the correlation matrix Amn- 
The matrix Amn is Hermitian, so it can be diagonalized by 
a unitary transformation. 



Gn 



m,?i— 



A 



(22) 



where = '^m "^pmCim- In this case, the density matrix of 
the block must also verify 



Gmn = Tr(.gl„5„pL) = Vm^mn , (23) 

which implies that p^ is uncorrelated and it can be written as 

PL = Qi® ■ ■ ■ ® QL , (24) 
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where is the density matrix corresponding to the m-th 
fermionic mode excited by g|„. 

In order to determine the eigenvalues of the density matrix 



of one mode, let us express the g„ 
its matrix representation. That is. 



gj^ and gm operators in 



and 



9rn = 





1 



I - Oir. 



1 





(25) 



(26) 



where a,„ and /?„, are the matrix elements of g,n that we want 
to determine. It is easy to see that = 0, since 

{g,n) = Tr(,9„iPL) = /3,„ = . (27) 

Moreover, rewriting Eq. |23]in terms of these matrices 



TT^idLgmPL) = Tr 



1 W a™ 
j I 1 - a„ 



(28) 

we realize that am ~ Vm. 

The entanglement entropy between the block and the rest 
of the system is therefore. 



(29) 



1=1 



where H2{x) = —x\ogx — (1 — a;)log(l — x) denotes the 
binary entropy. 

Summing up, the three steps that we have to follow in order 
to compute the entanglement entropy of the GS of a block of 
L spins for the XX model are: (;) to determine the correlation 
matrix A„in by evaluating Eq. ( fT9] l for L > m^n > 0, (;;) to 
diagonalize this correlation matrix and, with its eigenvalues, 
(///) to compute the entanglement entropy according to Eq. 



Let us emphasize that this method is computationally effi- 
cient, since its computational cost scales polynomially with 
the number of spins of the block 0{L^), whereas the Hilbert 
space of the problem has dimension 2^^. 

It is also necessary to recall a quite subtle point that we have 
skipped along our previous discussion. It turns out that there 
is no need to perform a final transformation back to spins, that 
is, there is no need to invert the Jordan- Wigner transformation. 
This is due to the fact that the coefficients of a given state are 
identical when written in terms of the spin basis or in terms of 
the fermionic a; operators. More precisely, 

H')^ J2 C''''-''"\ii,^2,...,iu) (30) 

il ,'12 : ■ ■■■in 

^ J2 C'-'--.^"(atrH4p...(at)'"|vac). 

Thus, the same coefficients appear in the ket, either when writ- 
ten in the initial spin basis, or when expressed as creation op- 
erators in the fermionic vacuum, |vac). Consequently, the re- 
duced density matrix entropy of entanglement is identical for 
both expressions. 



Finally, let us mention that the computation of the geomet- 
ric entropy of Gaussian systems have been systematized in 
Refs. iflql . In particular, it is shown that for solvable fermionic 
and bosonic lattice systems, the reduced density matrices can 
be determined from the properties of the correlation functions. 
This subject is reviewed in Ref. iITtIi in this special issue. 



E. Scaling of the entropy 

It is now easy to compute the entanglement entropy of the 
ground state of the XX model for arbitrary values of the block 
size L and magnetic field A. 




250 



FIG. 2: Entropy of the reduced density matrix of L spins for the XX 
model in the limit A'^ — > oo, for two different values of the external 
magnetic field A. The maximum entropy is reached when there is 
no applied external field (A = 0). The entropy decreases while the 
magnetic field increases until A = 2 when the system reaches the 
ferromagnetic limit and the ground state is a product state in the spin 
basis. 

In Fig. [21 we show how the entropy of the reduced density 
matrix of a block of L spins varies with L for different values 
of the magnetic field. The maximum entropy is reached for 
A = 0. In particular, we recover the result in Ref. ^ and see 
that for A = the leading scaling of the entropy is perfectly 
fitted by a logarithm. 



^L = ^ logs L 



(31) 



where a is a constant that was determined analytically in Ref. 

As we increase the magnetic field, but it is still less than 
2, the entropy decreases although it keeps its logarithmic be- 
haviour with L. When A > 2, the entropy saturates to zero, 
since the ground state is already in a product state in the spin 
basis corresponding to the ferromagnetic phase, J^^ 

The relation between logarithmic scaling and entropy is 
confirmed by similar computations in different models. The 
general result is that entanglement entropy obeys a logarith- 
mic scaling relation at critical points, that is when the system 
is at a phase transition, whereas a saturation of entanglement 
is found away from criticality. This universal logarithmic at 
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critical points must emerge from the basic symmetry that char- 
acterizes phase transitions, namely conformal invariance. We 
shall come to this developments in the next section. 

Let us mention that the scaling of entanglement entropy 
was formerly studied in the context of quantum field theory 
and black-hole physics. There, the system sits in higher di- 
mensions. The entanglement entropy scales following an area 
law that we shall discuss later on. Yet, it is important to note 
that one-dimensional systems are an exception to the area law. 
Entanglement pervades the system at any distance, not staying 
just at the point-like borders of a block. 

Summing up, we have seen that the scaling entropy is a 
good witness for quantum phase transitions. Many other stud- 
ies of different measures of entanglement at quantum phase 
transitions have been presented recently. Let us here mention 
that in Refs. lfl9[ I20I1 . quantum phase transitions are charac- 
terized in terms of the overlap (fideUty) function between two 
ground states obtained for two close values of external param- 
eters. When crossing the critical point a peak of the fidelity is 
observed. 



F. Entanglement entropy and Toeplitz determinant 

Before finishing this section, we would like to sketch how 
the particular structure of the correlation matrix in Eq. ( fT9] l 
allows us to derive an analytical expression for the scaling law 
of the entanglement entropy. This result is presented in Ref 
El. 

In order to obtain an analytical expression for the entangle- 
ment entropy, let us introduce the function 



(32) 



where A{fi) = ^II — A, II is the identity matrix of dimension 
L, and A is the correlation matrix defined in Eq. ( fT9] l. If we 
express the matrix A in its diagonal form, we trivially have 



(33) 



m — l 



where u„i, with m = 1, . . . , L, are the eigenvalues of A. 
Then, according to Cauchy residue theorem, the entanglement 
entropy Sl can be expressed in terms of an integral in the fi- 
complex plane as follows 



K <f e(l + e, fi)d In ^^(m) 
0+ s~*o+ 27ri /^(, ,5) 



Sl = lim lim - — : 

e 

L 



(34) 



m— 1 



where c(e, S) is a closed path that encircles all zeros of (/i) 
and e(l + e, ^) is an arbitrary function that is analytic in the 
contour c(e, S) and verifies e(l, i^m) = H2{vm) V m. 

Thus, if we could obtain an analytical expression for the 
Dl{^) function, we would be able to get a closed analytical 
result for the entanglement entropy. 



Notice that both the correlation matrix Amn, defined in Eq. 
( fT9] l, and A are Toeplitz matrices, that is to say, matrices in 
which each descending diagonal from left to right is constant. 



A = 



( /o /-I 
/i /o 



L-1 



/o / 



(35) 



1 sin kcTn 



where, in this case, /„ 

The asymptotic behaviour (when L — > 00) of the deter- 
minant of Toeplitz matrices has been widely studied in many 
cases, giving the famous Fisher-Hartwig conjecture (see Refs. 
Ill 23, 24, 25J). In our particular case, the expression for 
the determinant of A was proven in Ref. ll23ll and, therefore, it 
is a theorem instead of a conjecture. In this way, we may insert 
this result in Eq. ( [34l i, perform the corresponding complex in- 
tegral and obtain the asymptotic analytical expression for the 
entanglement entropy of the XX model. This computation is 
presented explicitly in Ref. [18] with the final result. 



where 



1 



■InL 



■ In 1 




(36) 



cosh(t/2) 1 



tsinli (t/2) 2sinh- 



'(V2)J 



(37) 

Indeed, we realize that this analytical expression for the scal- 
ing of the entanglement entropy is compatible with the nu- 
merical fit of Eq. ( [3TT i and, moreover, it fixes the value of the 
independent term. 

This procedure is also used to obtain an analytical expres- 
sion for the entanglement entropy of the XY model in Refs. 
ll26ll . In Ref. ll27ll . the scaling of the Renyi entropy is deter- 
mined for the XY model in terms of Klein's elliptic A - func- 
tion showing a perfect agreement with the previous results in 
the particular case in which the Renyi entropy becomes the 
von Neumann entropy. 



III. SCALING OF ENTANGLEMENT 

The logarithmic scaling law that entanglement entropy 
obeys in the critical regime is a sign of the conformal sym- 
metry of the system. For second order phase transitions, the 
correlation length diverges and the system becomes scale in- 
variant. This scaling symmetry gets enlarged to the conformal 
group li28ll which, in the case of on-dimensional systems, al- 
lows for a very precise characterization of the operator struc- 
ture of the underlying theory. The development of confor- 
mal field theory is a remarkable achievement that we cannot 
present in this short review ll29il30ll . 
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A. One-dimensional systems 

For a conformal theory in 1+1 dimensions, the scaling be- 
haviour of the entropy was proven to be logarithmic in Ref. 
lIsTll . The general result reads 

Sl ~ ^ log2 L , (38) 

where c and c are the so called central charges for the holo- 
morphic and anti-holomorphic sectors of the conformal field 
theory. These central charges classify conformal field theo- 
ries and are universal quantities which depend only on basic 
properties of the system, like effective degrees of freedom of 
the theory, symmetries or spatial dimensions, and they are in- 
dependent of the microscopic details of the model. For free 
bosons c = 1, whereas for free fermions c ~ 1/2. 

This result matches perfectly our geometric entropy com- 
putation of the critical XX model. In this case, the central 
charge c = c = 1 and the model is seen to belong to the free 
boson universality class. 

The previous result of Eq. dSSl l was further elaborated and 
extended to finite s ystems, finite temperature and disjoint re- 
gions in Ref. ll3Tll32l [3311 . For instance, the scaling of the 
entropy for a system with periodic boundary conditions reads 

^A-^logf— sin^) +4, (39) 
3 \na L J 

whereas for the open boundary conditions case is 

5^ = ^logf-sin^)+5;. (40) 
6 \na L J 

In Ref. ll34ll . the scaling of the entropy of a conformal semi- 
infinite chain is presented. In Ref. [135:1, conformal symmetry 
is further exploited and an analytical computation of the distri- 
bution of eigenvalues of the reduced density matrix of a block 
in a one-dimensional systems is presented. 

Let us finally mention that the scaling of entanglement have 
been also studied for other entanglement measures by means 
of conformal field theory. In particular, in Ref. ll36l[37ll . it is 
shown that the single copy entanglement scales as 

2 

E^{pL)^'-\ogL~j,-^ + 0{l/L). (41) 
6 6 log L 

Note that entropy sub-leading corrections to scaling are sup- 
pressed as 1/L whereas single copy entanglement suffers 
from 1/logL modifications. This makes the numerical ap- 
proach to the latter more difficult. 

B. Conformal field theory and central charge 

We stated above that the central charge is a quantity that 
characterizes the universality class of a quantum phase tran- 
sition. We also mentioned in the introduction that a possible 
figure of merit for entanglement could be constructed from 
correlation functions made out some operator which is always 



present in any theory. Let us now see that both ideas merge 
naturally. 

In 1+1 dimensions, conformal field theories are classified 
by the representations of the conformal group fl^. The oper- 
ators of the theory fall into a structure of highest weight oper- 
ators and its descendants. Each highest weight operator car- 
ries some specific scaling dimensions which dictates those of 
its descendants. The operators close an algebra implemented 
into the operator product expansion. One operator is particu- 
larly important: the energy-momentum tensor T^^, which is a 
descendant of the identity. It is convenient to introduce holo- 
morphic and anti-holomorphic indices defined by the combi- 
nations T = Tzz and T = Tzz where 2 = + ix^ and 
z ~ — ix^. Denoting by |0) the vacuum state, the central 
charge of a conformal field theory is associated to the coeffi- 
cient of the correlator 

(0|T(2)T(0)|0) = ^, (42) 

and the analogous result for c in terms of the correlator 
(0|T'(z)T'(0)|0). A conformal field theory is characterized 
by its central charge, the scaling dimensions and the coeffi- 
cients of the operator product expansion. Furthermore, uni- 
tary theories with c < 1 only exist for discrete values of c 
and are called minimal models. The lowest lying theory cor- 
responds to c = ^ and represents the universality class of a 
free fermion. 

The central charge plays many roles in a conformal field 
theory. It was introduced above as the coefficient of a corre- 
lator of energy-momentum tensors, which means that it is an 
observable. The central charge also characterizes the response 
of a theory to a modification of the background metric where 
it is defined. Specifically, the scale anomaly associated to the 
lack of scale invariance produced by a non-trivial background 
metric is 

(0|T^'|0)=-^i?, (43) 

where R is the curvature of the background metric. This 
anomaly can also be seen as the emergence of a non-local ef- 
fective action when the field theory modes are integrated out 
in a curved background. 

Therefore, the central charge which appears as the coeffi- 
cient of the entanglement entropy is naturally related to the 
stress tensor, which is the operator that is guaranteed to exist 
in any theory. It is also possible to derive a direct relation be- 
tween entropy and the trace of the stress tensor as shown in 
the original Ref. jJlll . 

C. Area law 

The conformal group does not constrain the structure of the 
Hilbert space in spatial dimensions higher than one as much 
as it does in one dimension. Actually, the conformal group no 
longer brings an infinite number of conserved charges (as it 
does in one dimension) but becomes a finite group. 

Nevertheless, a geometric argument establishes the scaling 
behaviour of entropy. The basic idea goes as follows. Let us 
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consider a volume of spins (or any local degrees of freedom) 
contained in a larger space. For theories with local interac- 
tions, it is expected that entanglement will be created between 
the degrees of freedom that lie outside and inside the surface 
that encloses the volume we are considering. It follows that 
the entropy should naturally scale as an area law even if the 
model displays a finite correlation length. 

These arguments were put forward in the study of entan- 
glement entropy in quantum field theory as a possible source 
for black-hole entropy iM HI ES- 

Furthermore, the rela- 
tion between the entropy and the effective action in a curved 
background was developed in Ref. ||4lll . Let us mention these 
results. For general quantum field theories the entropy is a di- 
vergent extensive quantity in more than one spatial dimension 
obeying an area law 



Sr 



Cl 



L 



d>l. 



(44) 



where L'^ is the size of the volume, e stands for an ultraviolet 
regulator and the coefficient ci counts components of the field 
which is considered. This coefficient is again found in the 
effective action on a gravitational field and, thus, in the trace 
anomaly as a divergent term. A form for the former can be 
found as 



eff 



ds- 



„d/2 



/Co 
V s 



CiR + C2FsF + C2GsG 



(45) 

where sg acts as a ultraviolet regulator, R is the curvature, F 
the Weyl tensor and G the Euler density. The main concep- 
tual result to be retained is that entropy measures a very basic 
counting of degrees of freedom. Note that previous efforts to 
make a general c-theorem are all base on C2f an C2g^ not on 
Cl . In one spatial dimension, the effective action has a unique 
structure proportional to the central charge. That is, the cen- 
tral charge takes over all manifestations of the trace anomaly, 
at variance with the separate roles that appear in higher di- 
mensions. 

It is worth mentioning that computations done in massive 
theories in any number of dimensions show that Sji m ^ 
0) — Sl (m = 0) comes out to be ultraviolet finite [42!||. Actu- 
ally, the ultraviolet cut-off cancels in the computation. This is 
precisely what is needed to make the RG flow meaningful in 
such a case. This comment hints at the non-trivial issue about 
observability of the entropy. The standard prejudice is that the 
leading area law coefficient is not observable since it comes 
divided by a necessary ultraviolet cut-off. Yet, if such a co- 
efficient is also responsible for finite corrections, the situation 
may not be as trivial. 

A review on methods to calculate the entanglement entropy 
for free fields and some particular examples in two, three and 
more dimensions are presented in Ref. 14311 in this special is- 
sue. Further explicit computations of area law scaling of en- 
tropy in spin and harmonic systems in higher dimensions can 
be found in Refs. H IH, 111 IH]. A quite remarkable re- 
sult found in Ref. ll48ll is the fact that certain fermionic sys- 
tems may develop logarithmic violations of the area law, while 
keeping local interactions. This is somehow analogous of the 



logarithmic scaling in one-dimension. The system is more 
correlated than what is expected from pure geometrical argu- 
ments. In this respect, the leading term in the scaling of the 
entropy for fermionic systems was computed analytically as- 
suming the Widom conjecture in Ref. 114911 . This result was 
checked numerically for two critical fermionic models in Ref. 
ll46ll finding a good agreement. 

For other steps into a description of systems with two spa- 
tial dimensions in the framework of conformal field theory see 
Refs. llsollsill . For a class of critical models in two spatial di- 
mensions (including the quantum dimer model), it is found 
that S{pi) = 2fs{L/a) + cg\og{L/a) + 0(1), where L is 
the length of the boundary area, fs is an area law coefficient 
that is interpreted as a boundary free energy, and g is a co- 
efficient that depends on the geometric properties of the par- 
tition. That is, in addition to a non-universal area law, one 
finds a universal logarithmically divergent correction. For a 
further discussion of steps towards a full theory of entangle- 
ment entropies in d + 1 -dimensional conformal field theories, 
see Refs. SH. 

A particularly interesting issue is the holographic entangle- 
ment entropy that emerges from the AdS (anti-de-Sitter)/CFT 
correspondence. The AdS/CFT correspondence is the conjec- 
tured equivalence between a quantum gravity theory defined 
on one space, and a quantum field theory without gravity de- 
fined on the conformal boundary of this space. The entan- 
glement entropy of a region of the boundary in the conformal 
field theory is then related with the degrees of freedom of part 
of the AdS space in the dual gravity side. In Refs. p U 15211 . 
this relation is established explicitly and, in Ref. ll53ll in this 
special issue, the recent progress on this topic is presented. 

Let us also mention the line of research that deals with topo- 
logical entropy. Some Hamiltonians produce states such that a 
combination of geometric entropies exactly cancels the domi- 
nant area law. Then, a topological entropy term characterizes 
the system iBill . This subject is nowadays a large field of re- 
search that we cannot include in the present review. In this 
respect, a review on the scaling of the entanglement entropy 
of 2D quantum systems in a state with topological order is 
presented in Ref. [55(1 in this special issue. 



rV. OTHER MODELS 

We can find in the literature the computation of the scaling 
of the entanglement entropy for other spin models. In XY and 
XXZ models, this logarithmic scaling will confirm the role 
of the underlying conformal symmetry. We shall also discuss 
that in disordered systems, although the conformal symmetry 
is lost for one particular realization of the disorder, we recover 
the logarithmic scaling of the entropy with a different central 
charge of the corresponding homogeneous model, if we do 
the average over all the possible realizations of the disorder 
We shall also study the scaling of entropy in systems where 
the notion of geometric distance is lost. This is the case of 
the Lipkin-Meshkov-Glick model, in which the logarithmic 
behaviour of the entropy will be due to the equilibrium of a 
competition between the long range interactions, that try to in- 
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crease the entanglement, and the symmetries of the problem, 
that force the ground state to belong to a reduced subspace 
of the Hilbert space. A different case are those systems com- 
posed of itinerant particles. In particular, we will present the 
scahng of entropy of the Laughlin wave function. 



A. The XY model 

The XY model is defined as the XX model in Eq. (|7]l, 
adding an extra parameter 7 that determines the degree of 
anisotropy of spin-spin interaction in the XY plane. Its Hamil- 
tonian reads 

(46) 

where, as in the previous section, I labels the N spins, (t[' 
(/i = X, y, z) are the Pauli matrices and A is the transverse 
magnetic field in the z direction. This notation will be also 
followed for the rest of models that are going to be presented. 

Notice that if 7 = we recover the XX model, whereas if 
7 = 1, it becomes the quantum Ising model with a transverse 
magnetic field, with Hamiltonian 

i?isi„g = -^^Kaf+i + Aaf). (47) 

The XY model was solved in detail in Ref. 101 ■ In order to 
do this, the previous works on spin chains required to solve the 
XY Hamiltonian were reviewed. In concrete, the XY model 
without magnetic field was solved exactly in Ref. lf56ll . the 
spectrum of the XY model with magnetic field was computed 
in Ref. iIstIi . the correlation function for this model was ob- 
tained in Ref. ifssll . and the entropy Sl was computed in Ref. 

Later, an extent analytical analysis of the entropy of XY 
spin chain was presented in Ref. ll26ll . In this work, in a sim- 
ilar way as we have seen previously for the XX model, an 
analytical expression for the scaling of the entanglement en- 
tropy is determined for the XY model by means of Toeplitz 
determinants. 

The XY model with 7 7^ is critical for A = 1. In this 
case, the entropy of a block scales as 

SxY{L) = \\og^L + a[-f), (48) 
6 

where 0(7) is a function that only depends on 7. This en- 
tanglement behaviour corresponds to the scaling dictated by a 
conformal theory, Eq. (l38T l, with central charge c = 1/2. The 
XY model, therefore, falls into the free fermion universality 
class. 

In the non-critical case, that is for A 7^ 1, the entanglement 
entropy saturates to a constant. 

Let us mention that an exact relationship between the en- 
tropies of the XY model and the XX model has been found 
recently ll59ll . Using this relation it is possible to translate 



known results between the two models and obtain, among oth- 
ers, the additive constant of the entropy of the critical homo- 
geneous quantum Ising chain and the effective central charge 
of the random XY chain. 

Finally, with respect to the particular case of the Ising 
model, in Ref. ll60ll . the computation of the leading correc- 
tion to the bipartite entanglement entropy at large sub-system 
size, in integrable quantum field theories with diagonal scat- 
tering matrices is presented. This result is used to compute the 
exact value of the saturation in the Ising model and showed it 
to be in good agreement with numerical results. This work is 
reviewed in detail in Ref. 1I6TI1 in this special issue. 

B. The XXZ model 

The XXZ model consists of a chain of N spins with near- 
est neighbour interactions and an external magnetic field. Its 
Hamiltonian is given by, 

Hxxz = E + <Tf af^i + Aafaf+J + Aaf) , 

(49) 

where A is a parameter that controls the anisotropy in the z 
direction. 

As it happened for the 7 parameter of the XY model, the 
A parameter of the XXZ model has two particular interesting 
values. If A = 0, we trivially recover the XX model, and if 
A = 1, the system becomes the XXX model that has a fully 
isotropic interaction 

Hxxx = (^K^f+i + + ^i^t+i] + M') ■ 

(50) 

The XXZ model can be solved analytically by means of 
the Bethe Ansatz technique ll62ll . Bethe Ansatz takes profit 
of the two symmetries of the system to find its eigenstates. 
The first symmetry is the rotational symmetry respect to the 
z axis. It implies that the z-component of the total spin, 
Sz ~ ^I'^^i'^i^ must be conserved and, therefore, the 
Hamiltonian must be diagonal in boxes of constant Sz- The 
other symmetry is the translational invariance, that allow us 
to diagonalize these boxes using a kind of generalized Fourier 
transform. Once the ground state is obtained, the correlation 
functions can be computed in terms of certain determinants 
(see Ref. ll63ll ). This model is qualitatively different from the 
XY, since it presents a point of non-analyticity of the ground 
state energy for finite systems. In the XY model, the level 
crossing between the ground state and the first excited state 
only occurs in the thermodynamic limit. In this case, instead, 
the terms of the Hamiltonian a'fa'f^^ + a^af^-^, erf af^^j^ and 
commute and are independent of A and A, which implies 
that there will be an actual level crossing. 

Both the isotropic case and the anisotropic one for A = 
are solved in Refs. ||2,|63l- The phases of the system are found 
to be: 

• In the XXX model, Eq. (|5Q] l, there are two limit be- 
haviours. On one hand, when |A| > 2 the system is 
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gapped and it is in a product state in which all spins 
point at the direction of the magnetic field (ferromag- 
netic phase). On the other hand, for A = the mag- 
netization is zero and the system is in a entangled state 
(anti-ferromagnetic case). In the interval between these 
two cases 2 > A > the system is gap-less and, there- 
fore, critical. 

• With respect to the anisotropic case with magnetic field 
equal to zero, the system shows a gap-less phase in the 
1 > A > — 1 interval. Outside this interval, there is 
a gap between the ground and the first excited states. 
These two phases are separated by two phase transitions 
in A = 1 and A = — 1. The first one is a Kosterlitz- 
Thouless phase transition, while the second one is of 
first order. 

The scaling of the entanglement entropy is presented in 
Refs. 10, lEsJ. These numerical results show that the entan- 
glement entropy behaviour converges to a logarithmic scaling 
as the size of the system increases, if the system is critical. On 
the contrary, if the system is not in a critical phase, the entropy 
saturates to a constant value. In particular, in the isotropic 
model without magnetic field, the entropy scales as 

SL^\\og^L, (51) 

which means that the XXX model with A = has central 
charge c = 1 and falls into the universality class of a free 
boson. 

Finally, let us mention that, recently, analytic expressions 
for reduced density matrices, several correlation functions and 
the entanglement entropy of small blocks (up to 6 spins) have 
been found for the XXZ model with A = 1/2 (see Refs. iH] 
and iH). 

C. Disordered models 

So far, we have only considered translational invariant sys- 
tems. This symmetry plus the scaling invariance at a critical 
point produces the conformal symmetry that implies universal 
properties of the scaling of entanglement. Nevertheless, natu- 
ral systems exhibit a certain amount of disorder due to impuri- 
ties and imperfections of the system. This disorder breaks the 
translational symmetry and we wonder what happens with the 
scaling of the entropy taking into account that the conformal 
invariance is lost. 

This question was addressed in Ref. ll68ll . They computed 
analytically the block entropy for the Heisenberg, XX and 
quantum Ising models with random nearest-neighbour cou- 
plings under the hypothesis of strong disorder by means of the 
real space renormalization group technique. This approach 
was introduced in Ref. ll69ll and was generalized in Ref. iFtoIi . 

This strong disorder hypothesis assumes that if one takes 
the strongest coupling of the chain, its neighbours are much 
weaker than it. Thus, it is possible to diagonalize this 
strongest bond independently of the rest of the system, project 
the system onto the ground state of this subspace (a singlet for 



the previous models) and then perform perturbation theory re- 
spect to the neighbour couplings. The final result is that two 
sites have been eliminated and the Hamiltonian energy scale 
has been reduced. This process can be iterated until we arrive 
at the ground state of the system which is a random singlet 
phase, that is to say, a set of singlets connected randomly and 
for arbitrarily long distances. 

Notice that although this method is not correct when ap- 
plied to a system with weak disorder, it becomes asymptoti- 
cally correct at large distances ll69ll . 

For a particular realization of the disorder, the translational 
symmetry of the system is broken and, therefore, the confor- 
mal symmetry too. Hence, the scaling of the entanglement 
entropy of this realization of the disorder will not be logarith- 
mic, but fluctuating. 

In Ref. ll68ll . it was shown that although the conformal sym- 
metry is broken, if we take the average over all the realizations 
of the disorder the entropy keeps scaling logarithmically with 
an effective central charged c = cln 2, where c is the central 
charge for the same model but without disorder. This result 
has been further checked numerically both for the XX model 
in Ref. 17 111 and for the Heisenberg model in Ref. ll72ll . 

In Ref. ItiIi . the disordered XX spin-i chain with periodic 
boundary conditions and positive random spin couplings cho- 
sen in a flat uniform distribution within the interval [0, 1] was 
studied. The magnetic field was set to zero. It was shown 
that for a block large enough (larger than 20 spins), the en- 
tropy scales logarithmically according to ll68ll . using around 
Iff^ samples for N = 500, 1000, 2000 and 2 x lO"* samples 
for 100 < N < 400, in order to do the average over all the 
possible realizations of the disorder 

The same result was shown for the Heisenberg model in 
Ref. ll72ll . In this work, a uniform distribution in the inter- 
val [0, 1] for the couplings between the spins was also chosen. 
For a system of = 50 and after averaging the entanglement 
entropy over lO''^ different configurations of disorder, the loga- 
rithmic scaling of the entropy with an effective central charge 
c = cln 2 is recovered. Let us point out that these one di- 
mensional systems are particular cases of chains of quantum 
group (or q-deformed) spins studied in Ref. iItsIi . It is also in- 
teresting to mention that this robustness of the entanglement 
scaling respect to the disorder is not kept for other models like 
the Bose-Hubbard model (see Ref. H). 

In the case of higher dimensions, the scaling of the entan- 
glement entropy in the 2D random Ising model was studied 
in Refs. ll75ll76fl . In particular, in Ref. ifTql . the entanglement 
entropy of a L x L region located in the centre of a square 
lattice which is governed by the Hamiltonian 

^ = -E^'^<^I-E^^<' (52) 

was computed. The Ising couplings Jy and the transverse 
magnetic fields Aj take random values given by the uniform 
probability distributions in the intervals [0, 1] and [0, Aq] re- 
spectively. By means of a generalized version for 2 dimen- 
sions of the real space renormalization group, it was found 
that the critical field is at AS = 5.37 ± 0.03, and for both criti- 
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cal and non-critical Ao the entropy scaling fulfils the area law: 
S{L) ^ L in the leading term. 

Let us mention some disordered spin systems have also 
been studied from the fidelity susceptibility point of view in 
Ref. IttIi . Finally, it is interesting to point out that, in other 
systems, the translational invariance is not broken by means 
of random couplings but due to a quantum impurity or a phys- 
ical boundary. The behaviour of the entanglement entropy in 
this kind of systems is reviewed in Ref. IITSll in this special 
issue. 



D. The Lipkin-Meshkov-Glick model 

The Lipkin-Meshkov-Glick model was proposed in Ref. 
ll79ll . Unlike the previous models we have considered, where 
the spins had short range interactions, in the LMG model, each 
spin interacts with all the spins of the system with the same 
coupling strength. This system is described by the Hamilto- 
nian 

Hlmg = E + ^^M) - E (53) 

i<j i 

Notice the apparent similarity between this model and the XY 
model in Eq. ( l46b . The essential difference between them is 
that while in the XY Hamiltonian the interaction only takes 
place between nearest neighbours, in the LMG model, all 
spins interact among themselves. This highly symmetric inter- 
action pattern forces the loss of the notion of geometry, since 
there is no distance between the spins. This implies that it 
no longer makes sense to define a block of L spins as a set 
of L contiguous spins or to study decays of the correlations 
between two spins. 

As in previous cases, our aim is to study the scaling proper- 
ties of the entanglement entropy for the ground state reduced 
density matrix of a block of L spins respect to the rest of Nl 
spins. We face a somewhat contradictory situation. On one 
hand, we expect that the non-local connectivity of the inter- 
actions would produced a ground state more entangled than 
those that emerge from nearest neighbour interaction models. 
On the other hand, the symmetry of the Hamiltonian implies 
that all the spins must be indistinguishable in the ground state, 
therefore, it must belong to a symmetric subspace, which re- 
stricts its entanglement. The explicit computation will clarify 
this issue. 

The Hamiltonian ( |53] ) can be written in terms of the total 
spin operators Sa = J2i as 

i7 = -1(1 + 7) (S^ - 5? - iV/2) - 2XSz 

-^(1-7) (5! + 5^) , (54) 

where S± are the ladder angular momentum operators. In Eq. 
( l54l i. we realize that [S^, H] = and, therefore, we can di- 
agonalize the Hamiltonian in boxes of constant S. From Eq. 
(|54] i, it is easy to see that the ground state must belong to the 
subspace of S* = N/2. Then, we have to span this subspace 



in terms of a basis \N/2, M) fully symmetric under the per- 
mutation group and eigenstates of and Sz- These states 
I M) are called Dicke states. 

Notice that the restricted subspace where the ground state 
must live due to the symmetries of the Hamiltonian will 
limit the scaling of the entanglement entropy of a block of 
L spins with respect to the remaining N^. As the ground 
state reduced density matrix is spanned by the set of (L + 1) 
Dicke states, the entropy of entanglement obeys the constrain 
Sl,n < log2(^ + 1) for all ^ and N, where the upper 
bound corresponds to the entropy of the maximally mixed 
state pl,n = ^/{L + 1) in the Dicke basis. This argument 
implies that entanglement, as measured by the Von Neumann 
entropy, cannot grow faster than the typical logarithmic scal- 
ing observed in the previous cases. 

Both the ground state and the entanglement entropy were 
computed for the LMG model in Ref. llsoll . For the isotropic 
case (7 = 1) and in the thermodynamic limit {N, L ^ 1), 
Hlmg is diagonal in the Dicke basis. Then, for A > 1, the 
entanglement entropy is strictly zero since the ground state is 
in a fully polarized product state. Instead, if 1 > A > 0, we 
recover the logarithmic scaling of the entropy, 

SL,N{\l^^)-\^og[L{N-L)/N]. (55) 

Although the kind of scaling does not depend on the strength 
of the magnetic field, its absolute value is smaller for weaker 
magnetic fields, according to the equation 

Sl.n{\i^ 1) - Sl.n{\ = 0,7 = 1) ~ ^log (1 - A^), 

_ (56) 

and thus diverges, at fixed L and N, in the limit h ^ 1^. 

In the anisotropic case, we can study the limits of very 
strong and very weak magnetic fields. On one hand, when 
A ^ 00, the GS is in the product state and therefore 

is not entangled. In the thermodynamic limit, this state is also 
the ground state just for A > 1. On the other hand, for A ^ 
the entanglement entropy saturates and goes to a constant that 
depends on 7. In the particular case of 7 = 0, the ground state 
is degenerate and lives in the subspace generated by the states 
rii and rii where |^) and |^) are the eigenstates 
of the cr^' operator In practice, this degeneration would be 
broken by any perturbation of the environment. 

These two different phases suggests the existence of a quan- 
tum phase transition between A ^ 1 and A <C 1. In partic- 
ular, it has, numerically, been checked in Ref. ifsoll that, in 
the thermodynamic limit, the entanglement entropy displays a 
logarithmic divergence around Ac = 1 according to the law 

5L,iv(A,7)^-log|l-A|. (57) 

Indeed, it is shown that at A = 1 the entropy scales logarith- 
mically with a coefficient that depends on 7. However, in the 
thermodynamic limit, this coefficient is independent of 7 and 
takes a value closed to 1/3. In Ref. Isill . the previous relation, 
Eq.|57] is computed analytically obtaining the same result and 
fixing the coefficient to 1/3. In this same work, the finite size 
corrections to the scaling of entropy are also studied. 
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Although the behaviour of entanglement is very similar to 
the XY model, that is to say, it scales logarithmically in the 
critical point and saturates to a constant in the non-critical 
phase, the reasons of these scaling laws are different. In the 
XY model, entanglement is limited by the facts that interac- 
tions are local and the system is translationally invariant. At 
the critical point, the correlation length becomes infinite, the 
system is conformal symmetric, and the logarithmic scaling 
of the entanglement entropy appears as a manifestation of this 
symmetry. Instead, in the LMG model, the long range in- 
teractions should allow for larger correlations, that is, larger 
entanglement. Nevertheless, the symmetries of the system re- 
strict the subspace where the GS must belong and, therefore, 
the scaling law of entanglement. The final result is the same 
logarithmic scaling law but which has nothing to do with any 
underlying conformal symmetry. 

Finally, let us mention that other analytical calculations of 
the spectrum of the LMG model both in the thermodynamic 
limit and finite size case have appeared recently [82,1 . More- 
over, the entanglement entropy for general free bosonic two- 
mode models is presented in Ref. |[83i. In particular, a com- 
plete classification of the possible scaling behaviours for the 
entanglement entropy in the related collective models as the 
LMG, the Dicke model, or the Lieb-Mattis model is obtained. 



where Zj = Xj +iyj,j = 1, . . . ,n stands for the position of 
the j-th particle. It describes fractional quantum Hall state at 
a filling fraction i/ = 1/m, where m is an integer number 

In particular, in Ref. IstIi . the entanglement entropy of k 
particles respect the rest of the system is computed for the 
Laughlin wave function with filling fraction one 

Sn,k = l0g2 (^^ ■ (59) 

Notice that, in this case, although the state also belongs to 
a completely anti-symmetric subspace, the entanglement en- 
tropy of half a system grows linearly with the number of par- 
ticles. 

In Ref. Issll . these ideas are extended, and the particle en- 
tanglement, defined as the entanglement between two subsets 
of particles making up the system, is studied. The general 
structure of particle entanglement in many-fermion ground 
states, analogous to the area law for the more usually stud- 
ied entanglement between spatial regions, are also formulated, 
and the basic properties of particle entanglement are uncov- 
ered by considering relatively simple itinerant models. All 
these ideas are widely reviewed in Ref. |[89ll in this special 
issue. 



E. Particle entanglement 

In a similar way than LMG model, where the notion of 
distance was lost, one can try to compute the entanglement 
entropy in systems of moving fermions and bosons. In such 
itinerant systems, as the particles are indistinguishable, mov- 
ing and partially de-localized, it is not obvious to define the 
geometric entropy. 

What we, indeed, can compute is the von Neumann entropy 
for any subset of particles for a system of indistinguishable 
particles in the state ^(ri, . . . , r„). Notice that, in this case, 
this von Neumann entropy cannot be interpreted as the num- 
ber of distillable EPR pairs. Due to the symmetrization, it is 
impossible to associate a label with the particles and perform 
the appropriate distillation operations. This is a subtle differ- 
ence respect to the LMG model. 

A particular interesting physical system is the Fractional 
Quantum Hall Effect (FQHE) Q. Although a complete un- 
derstanding of it is still missing, it is commonly believed that 
the interactions between the particles are essentially respon- 
sible for the strange states of matter that the 2D electron gas 
shows at some particular values of the transverse magnetic 
field. These states would present a new kind of order called 
topological order and their quasi-particle excitations are nei- 
ther bosons nor fermions, but anyons, that is to say, quasi- 
particles with any-statistics Issll . In this respect, in 1983, 
Laughlin proposed an Ansatz for the wave function of the 
ground state of the system Js^. This wave function is defined 
by 

■■■,zn)- X{{z^ - z,re- 1 , (58) 

i<3 



V. RENORMALIZATION OF ENTANGLEMENT 

A natural question arising within the study of entanglement 
in quantum system is how entanglement evolves along Renor- 
malization Group (RG) trajectories. We shall now address this 
issue discussing first the RG of quantum states and, then, the 
study of particular systems. 



A. Renormalization of quantum states 

It is customary to present RG transformations on Hamil- 
tonians or observables. In general, a Hamiltonian is de- 
scribed by a set of coupling constants times operators H = 
J2i9^^i- This set of operators may be infinite, including 
relevant, marginal and irrelevant operators or, as in the case 
of renormalizable quantum field theories in the continuum, 
it may reduce to a finite set of relevant and marginal oper- 
ators. Then, upon coarse graining of short-distances and an 
adequate rescaling, the system is described by a new set of 
coupling constants. So to speak, the operator algebra acts as a 
basis. The concept of RG trajectory corresponds to analysing 
observables along the flow ^ ~ df^' where the beta func- 
tions correspond to (3i = ^ are related to the change of the 
coupling constants as the coarse graining proceeds. 

Yet, RG transformations can be understood as an action on 
any quantum state, regardless of its relation to any Hamilto- 
nian ll90ll . Coarse graining is independent of any dynamics. 
This RG procedure on quantum states is not presented as com- 
mon lore since explicit knowledge of e. g. the ground state of 
a system is not available in general. Let us address this issue. 



13 



The basic idea to perform RG on states is to produce a 
coarse graining of short-distance degrees of freedom, fol- 
lowed by a clever choice of local basis to retain the long- 
distance information which is retained in an optimal way. Let 
us take a quantum state ipo and determine its RG transformed 
, ipQ, as follows. We pairwise group the sites in the system and 
define a coarse-graining transformation for every pair of local 
d-dimensional basis states, e.g. for the sites 2j and 2j + 1, 
as \p)2j\q)2j+i = \PQ)j- This transformation yields -^o — ^ V'- 
Then we have ipQ = U ^ . . .(E)U\tl') , where the (P x d? unitary 
matrix U performs the change of representative in the coarse- 
grained space. Note that the matrix U is non-local as seen 
from the 2j and 2j + 1 sites. Some local information is now 
washed out, while preserving all the quantum correlations re- 
lating the coarse-grained block to other ones. 

Operators also get coarse-grained along the above trans- 
formation. Take for instance an operator acting on one local 
Hilbert space, e.g. 02j ■ Expectation values must remain un- 
changed, 

(V-olOsjIV^o) = (V-olO^-lV^) , (60) 

which leads to 

0[ = U{02j®l2Hi)U\ (61) 

where / is the identity matrix. To complete a RG transforma- 
tion we simply need to rescale distances, i.e., to double the 
lattice spacing. 

This analysis can be made completely explicit in the case 
of states which are described as a matrix product state ll90ll . 
There, the above transformation amount to a flow on the ma- 
trices that represent the state. In turn, a flow related to the 
transfer matrix can be computed. Explicit irreversibility of 
RG flows and the characterization of critical points followed 
from the flow on this transfer matrix. 



B. Irreversibility of RG flows 

We may as well return to the standard construction of RG 
transformations on Hamiltonians and perform a detailed study 
in some particular case. For instance, we may consider the 
quantum Ising model in a transverse field A. It is known that 
the parameter A provides a relevant deformation of the model, 
departing form its critical value A* = 1. For instance, the de- 
parture that makes A > 1 get larger and larger corresponds to 
the increase of the mass of the underlying fermionic descrip- 
tion. 

An analysis of this RG trajectory can be illustrated using 
Fig.[3]( see Refs. MM)- 

This result shows that RG tra- 
jectories are monotonically irreversible as dictated by the c- 
theorem in 1h-1 dimensions. Furthermore, it can be seen that 
the ground state obeys majorization relations. That is, irre- 
versibility is orchestrated at a very refined level, since the 
reshuffling of the ground state obeys an exponential set of or- 
dering relations ll44ll . 

Irreversibility for the entanglement entropy should then 
be obtained as a fundamental theorem, equivalent to the c- 
theorem which is usually formulated in terms of the stress 
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FIG. 3: Entropy of entanglement is shown to decrease monotonically 
along the RG trajectory that takes the external magnetic field A away 
from its critical value A* = 1. Towards the left the flow takes the 
system to a GHZ-line state whereas, towards the right, the system is 
a product state. 

tensor. This was indeed done in Ref. ll92ll . Once, the rela- 
tion between entropy and the properties of the stress tensor 
are made apparent. 

VI. DYNAMICS OF ENTANGLEMENT 

So far, we have studied the properties of entanglement en- 
tropy of the ground state of the system. Next, we would like to 
analyse how entanglement evolves in time when the system is 
prepared in a state that is not an eigenstate of the Hamiltonian. 

A. Time evolution of the block entanglement entropy 

In Ref. 1I93I1 . the time evolution of the entropy of entan- 
glement of a block of L spins in a one-dimensional system is 
studied. It is considered a system prepared in a pure state IV^o)- 
which corresponds to an eigenstate of H{\o) with Ao 7^ A. 
Then, for example, at time t ~ 0, the parameter is suddenly 
quenched from Aq to A. In general, {ipo) will not be an eigen- 
state of -ff(A), and thus the system will evolve according to 
the equations of motion given by H{X). In this work, two 
computations are performed; one based on conformal field 
theory and the other on a particular solvable spin model, the 
Ising model. In the first case, the path integral formulation 
and the CFT are used in order to calculate the time evolution 
of the entanglement entropy of a high energy state of the sys- 
tem which is not an eigenstate. Then, one has to assume that 
the Hamiltonian is critical in order to make the theory confor- 
mally invariant. Instead, in the Ising model case, it is possible 
to perform calculations starting from a variety of initial states, 
considering both critical and non-critical regimes. 

In both calculations, the entanglement entropy increases 
linearly with time t (after transients die away in the lattice 
case), up to t* = L/2, in units where the maximum propa- 
gation speed of excitations is taken to be unity. For t ^ t*, 
Shit) ^ L saturates at an asymptotic value. This behaviour 
can be summarized in the following equation: 

SL{t)^[l l^l^ . (62) 
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This behaviour of the entanglement entropy has been 
checked in several lattice models both analytically and nu- 
merically In particular, in Ref. ll94ll . the 
previous results are provided analytically using Toeplitz ma- 
trix representation and multidimensional phase methods for 
the XY model and considering large blocks. 

In Ref. 1I93I1, a simple interpretation of this behaviour is 
proposed in terms of quasi-particles excitations emitted from 
the initial state at t = and freely propagating with velocity 
V < 1. The idea is that at, t = and at many points of the 
chain, a pair of entangled quasi-particles begin to propagate 
in opposite directions at some constant velocity v that we will 
consider 1 for simplicity (see Fig|4]i. The entanglement be- 
tween the block of L spins and the rest of the system at an 
arbitrary time is given by the number of pairs that have one 
quasi-particle in the block while the other is outside. Thus, 
the entanglement entropy increases Unearly with time until it 
saturates when the excitations that started in the middle of the 
block arrive at its boundary. 

All the previous results are explained in detail in the Ref. 
Issll . where, apart from quantum quenches, a general con- 
formal field theory approach to entanglement entropy is re- 
viewed. 
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FIG. 4: Schematic representation of the dynamics of block entropy. 
Entangled particles are emitted from the region A, they will con- 
tribute to the block entropy as long as one of the two particles ends 
in the region B [from 19311 1. 

Let us point out that this increase of the entanglement en- 
tropy is unrelated to the second law of thermodynamics. En- 
tanglement entropy can decrease or even oscillate in standard 
time evolution. 

Let us also mention that in Ref. ||72ll the dynamics of entan- 
glement was analysed for disordered systems, i. e. when the 
couplings between the spins take random values. In particular, 
the XXZ model with the couplings between the spins follow- 
ing a uniform distribution in the interval [0, 1] was studied. It 
turns out that, in the presence of disorder, the entanglement 
entropy does not increase linearly but logarithmically. This 
logarithmic behaviour does not follow from an extension of 
the argument for the clean case assuming a diffusive propa- 
gation of the excitations, but it requires some kind of entan- 
glement localization. This behaviour is also observed in Ref. 
119811 where the propagation of information through the disor- 



dered XY model is studied. In particular, both classical and 
quantum correlations are exponentially suppressed outside of 
an effective light-cone whose radius grows at most logarith- 
mically with time. 

B. Bounds for time evolution of the block entropy 

All these results a re co mpatible with the rigorous bounds 
found in R efs. ||99[ llOOll by means of the Lieb-Robinson 
bound i lOll l and its generalizations presented in Refs. 1 102]. 

The Lieb-Robinson bound states that the operator norm of 
the commutator of two operators Oa and O b that act on dif- 
ferent regions A, i? of a spin network with local interactions, 
hij (t), and in different times verifies 

\\[OA{t),OBm\\ < cN,n^n\\OA\\\\OB\\e-^ , (63) 

where L is the distance between A and B (the number of 
edges in the shortest path connecting A and E), N„iin = 
min{|A|, \B\} is the number of spins in the smallest of A 
and B, while c, u, ^ > are constants depending only on 
g = max(j maxt (t)|| and the architecture of the 
spin lattice. 

Thus, the Lieb Robinson bound, Eq. ( |63] |. tells us that the 
norm of the commutator of two operators at different times is 
exponentially small outside the light-cone given by the veloc- 
ity V that we can understand like the speed of sound. Notice 
that, by dimensional analysis, this velocity must be propor- 
tional to the energy scale g. It is interesting to point out that 
this result is also valid for the case of fermions or local Hamil- 
tonians with exponentially decaying interactions. 

In Ref. ll99ll . it is sh own, using the Lieb-Robinson bound 
and its generalizations ll02ll . that correlations and informa- 
tion are propagated at a finite velocity in a spin network with 
nearest-neighbour interactions. This is a non-trivial result 
since in non-relativistic quantum mechanics there doesn't ex- 
ist the notion of a light-cone, i. e. local operations could be 
used, in principle, to send information at long distances in ar- 
bitrary short times. 

Moreover, it is quantified the entanglement entropy that can 
be generated per unit of time between a block of spins and the 
rest of the system. In particular, it is found that 

SLit)-SLiO)<c*gPt (64) 

where c* ~ 1.9 is a constant and P is the perimeter of the 
block. Finally, let us m ention that all these results are com- 
plemented in Ref. ifiooli . 

C. Long range interactions 

The Lieb-Robinson bound is only valid for short range in- 
teractions. Then, it is interesting to study how does entan- 
glement evolve in systems wi th lo ng range interactions. This 
question is addressed in Ref. Il03ll . 

In general, systems with long range interactions are nu- 
merically intractable since, in them, the entanglement entropy 
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scales with the volume Sl ~ L. Nevertheless, in Ref. lll03ll . 
the interactions are restricted to Ising-type which allows to 
study both the static and the dynamical entanglement proper- 
ties of the system. 

It is considered a lattice composed by N spins that interact 
according to the Hamiltonian, 

^ = E J(l - ^i'^) ® (1 - ^i'') ' (65) 

k<l 

where the coefficients f{k, I), that describe the strength of the 
interaction between the spins I and k, obey a distance law, that 

istosay,/(fc,0 = /(|| fc-MD- 

It is assumed that the initial state is a product state of all 
spins pointing to the x-direction |4'o) = | -^)'^^ . In order to 
perform the time evolution of this state, a description in terms 
of Valence Bond SoHds (VBS) is used (see Ref. ifioill ). With 
this method, it is possible to calculate the reduced density op- 
erator of few particles for large systems (the computational 
time grows linearly with the whole size of the system but ex- 
ponentially with the size of the block). 

In concrete, it is studied for some fixed time t the scaling 
properties of entanglement of a system with algebraically de- 
caying interactions /(fc,Z)=|jfc — It turns out that for 
<x < 1/2 (strong long-range interactions) the entanglement 
grows unbounded and the correlations do not practically de- 
cay, while for a > 1 the system contains a bounded amount 
of entanglement and the correlations decay algebraically. 

The dynamics of entanglement are also studied. In the limit 
of an infinite chain, the entanglement entropy of any block 
saturates for large times (t oo) to its maximal value Sl = 
L in a similar way as in Eq. ( |62] |. 



VII. ENTANGLEMENT ALONG QUANTUM 
COMPUTATION 

It is known that slightly entangled quantur n sys t ems can be 
simulated efficiently in a classical computer lll05lll06ll . This 
implies that any quantum algorithm that would exponentially 
accelerate a classical computation must create, at some point, 
a highly entangled state. Otherwise, the quantum algorithm 
could be simulated efficiently in a classical computer. 

Next, we want to briefly study how the entanglement 
evolves along a computation. In order to do this, we will con- 
sider the three most common paradigms of quantum compu- 
tation: quantum circuits, adiabatic quantum computation, and 
one way quantum computing. 

A. Quantum circuits 

A quantum circuit is a sequence of unitary transf ormat ions 
(quantum gates) on a register of qubits (see Ref. lll07ll for 
a pedagogical introduction). An efficient quantum circuit is 
characterized by the fact that the number of elementary gates 
that form it only scales polynomially in the number of qubits 
of the register 



The study of en tangl emen t alon g a quantum circuit was ad- 
dressed in Refs. lIlOSll and lll09[l by means of majorization 
theory. In these works the introduction of entanglement in 
Shor's algorithm and the Grover's algorithm were analysed 
respectively. 

Let us remind the concept of Majorization relations, which 
is a more refined measure of ordering of probability distribu- 
tions than the usual entropy one. We say that a probability dis- 
tribution {pi} majorizes another probability distribution {qi} 
(written as p ^ q) if, and only if, 

k k 

E^'^^E* k = l...d-l, (66) 

1=1 1=1 

where d — 2^ is the number of possible outcomes and it will 
correspond to the dimension of the Hilbert space. 

These Majorization relations can be related to quantum cir- 
cuits in the following way: let \ipm) be the pure state rep- 
resenting the register in a quantum computer in the com- 
putational basis at an operating stage labeled by m = 
0,1 ... M — 1, where M is the total number of steps of 
the algorithm. We can naturally associate a set of sorted 
probabilities pi™' corresponding to the square modulus of 
the coefficients of the state in the computational basis (x E 
{|0...0),|0...01),...,|l...l)}). A quantum algorithm 
will be said to majorize step by step this probability distri- 
bution if 

^rn) ^ Vto = 1, . . . , Af . (67) 

In such a case, there will be a neat flow of probability di- 
rected to the values of highest weight, in a way that the proba- 
bility distribution will be steeper and steeper as the algorithm 
goes ahead. This implies that the state is becoming less en- 
tangled along the computation. Notice that the majorization 
relations are stricter than an inequality in the entanglement 
entropy, in such a way that the reverse statement is not true. 

In Ref. ll lOll . the step-by-step majorization was found in 
the known instances of fast and efficient algorithms, namely in 
the quantum Fourier transform, in Grover's algorithm, in the 
hidden affine function problem, in searching by quantum adi- 
abatic evolution and in deterministic quantum walks in con- 
tinuous time solving a classically hard problem. On the other 
hand, the optimal quantum algorithm for parity determination, 
which does not provide any computational speed-up, does not 
show step-by-step majorization. 

Recently, a new class of quantum algorithms have been 
presented. Those are exact circuit that faithfully reproduce 
the d ynam ics of strongly correlated many-body system. In 
Ref. Ill ill , the underlying quantum circuit that reproduces the 
physics of the XY Hamiltonian for N spins was obtained. The 
philosophy inspiring that circuit was to follow the steps of the 
analytical solution of that integrable model. Looking at the 
architecture of the circuit in Fig.|5] it is easy to realize that the 
entanglement between the two sets of contiguous N/2 spins 
is transmitted through the N/2 SWAP gates. Therefore, the 
maximum entanglement entropy between these two half's of 
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FIG. 5: Structure of the quantum circuit performing the exact diago- 
nalization of the XY Hamiltonian for 8 sites. The circuit follows the 
structure of a Bogoliubov transformation followed by a fast Fourier 
transform. Three types of gates are involved: type-B (responsible for 
the Bogoliubov transformation and depending on the external mag- 
netic field A and the anisotropy parameter 7), type-fSWAP (depicted 
as crosses and necessary to implement the anti-commuting proper- 
ties of fermions) and type-F (gates associated to the fast Fourier 
transform). Some initial gates have been eliminat ed si nce they only 
amount to some reordering of initial qubits [from 111 ih l. 

the system that this proposal may allow is N/2. This is be- 
cause the maximum entanglement that can generate a quan- 
tum gate that acts on two qubits is 1, that is, from a product 
state to a maximally entangled state (Bell basis). Thus, the 
scaling law of the entanglement entropy that this proposal will 
allow will be 

S{N/2) < N/2 . (68) 

Notice that, as we have seen in the previous sections, the en- 
tanglement entropy of the ground state in the XY model scales 
only logarithmically. The above circuit, then, can create much 
more entropy than what is present in the ground state. Yet, 
we have also discussed the fact that time evolution does cre- 
ate maximum entanglement. This, indeed, is what the above 
circuit achieves. This shows that the previous proposal is op- 
timal since it carries the minimum possible number of gates 
such that maximum entanglement can be created. 

Let u s als o add a final example on exact quantum circuits. 
In Ref. ll I2II . a quantum circuit that creates the Laughlin state 
(Eq. |58] | for an arbitrary number of particles (qudits) n in the 
case of filling fraction one is presented. The way in which 
entanglement grows along the circuit is also related to the 
amount of entanglement that each gate of the circuit can gen- 
erate. In the case of this Laughlin wave function, the depth of 
the circuit grows linearly with the number of qudits, so linear 
entanglement S ^ n can be supported by the circuit. This 
is precisely the entanglement that the Laughlin wave function 
with filling fraction one requires as shown in Ref. IstIi . 

The exact circuits we have discussed (XY and Laughlin 
states) are both able to create linear entanglement entropy. It 
is then impossible they can be simulated classically in an effi- 
cient way. 

B. Adiabatic quantum computation 



ing: 

1 . A quantum register is initially prepared on the ground 
state of a known initial Hamiltonian Hq. 

2. The system is then made to evolve adiabatically from 
this Hamiltonian to a new one Hp whose ground state 
codifies the solution to an e.g. NP-complete problem 

H{.s{t)) = (1 - s{t))Ho + s{t)Hp . (69) 

3. Slow evolution from s(t = 0) = to s(t = T) = 1 
guarantees that the system will not jump from the in- 
stantaneous ground state of the system to the first ex- 
cited state. 

Quantum adiabatic computation is proven efficient provided 
that the minimum gap along the adiabatic evolution is only 
polynomially small in the number of qubits. If this was not 
the case, the adiabatic computation would require an exponen- 
tially large time as measure in terms of the number of qubits 
in the register. 

Thus, according to the previous arguments, at some point 
of the adiabatic evolution of a hard quantum computation the 
system must be highly entangled, in a similar way as it hap- 
pened in the previous sections at the quantum phase transi- 
tions. This makes us expect some sort of quantum phase tran- 
sition for a concrete value Sc of the Hamiltonian, point that 
would be c harac terized by a minimum energy gap. 

In Ref. ill 14ll . adiabatic quantum computation is used to 
solve the NP-Complete Exact Cover problem that is a particu- 
lar case of the 3-SAT problem. It is defined as follows: given 
the n Boolean variables {xi}i=i^,,,n, Xi = 0, 1 V i, where i is 
the bit index, we define a clause C involving the three bits i, j 
and k by the constraint Xi + Xj + Xk = 1. There are only three 
assignments of the set of variables {.t^, Xj, Xfc} that satisfy 
this equation, namely, {1, 0, 0}, {0, 1, 0} and {0, 0, 1}. An /«- 
stance of the Exact Cover problem is a collection of clauses 
which involves different groups of three qubits. The problem 
is to find a string of bits {xi, X2 ■ ■ ■ ,Xn} which satisfies all 
the clauses. 

This problem can be mapped to findin g the ground state of 
a Hamiltonian Hp in the following way |l 14!]: given a clause 
C define the Hamiltonian associated to this clause as 

i/c = ^((l + <)(l + 'T|)(l + a|) (70) 

+ (1 - )(1 ~ a|)(l - al) + (1 - <Tf)(l - + 

+ (1 - <7f )(1 + a|)(l - al) + (1 + )(1 - a|)(l - , 

where cr^ |0) = |0), cr^ll) = — The quantum states of 
the computational basis that are eigenstates of He with zero 
eigenvalue (ground states) are the ones that correspond to the 
bit string which satisfies C, whereas the rest of the com- 
putational states are penalized with an energy equal to one. 
The problem Hamiltonian is constructed as the sum of all the 
Hamiltonians corresponding to all the clauses in the instance. 



The framework of adiabatic quantum computation (AQC) 
was introduced in Ref. 1 1 1 3ll . The idea of AQC is the follow- 



Hp — ^ He 

C tz instance 



(71) 
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FIG. 6: Average over 300 instances of the entanglement entropy be- 
tween two blocks of size 7i/2 as a function of the parameter s control- 
ling the adiabatic evolution. A peak appears for Sc ~ 0.7. The plot 
also shows the increase of the peak as the number of qubits grows 
n = 10, 12, 14 [from jH]. 



The ground state of this Hamiltonian corresponds to the quan- 
tum state whose bit string satisfies all the clauses. 

It is known that Exact Cover is a NP-complete problem, so 
it cannot be sol ved in a polynomial number of steps in a clas- 
sical computer ll 15ll . This makes the Exact Cover problem, 
particularly interesting, since if we had an algorithm to effi- 
ciently solve Exact Cover, we co uld also solve all problems in 
the much l arger NP family ill 16ll . 

In Ref. ill 17ll . the evolution of the entanglement properties 
of the system are studied in order to see the expected sign 
of a quantum phase transition. 300 random instances for the 
Exact Cover are generated with only one possible satisfying 
assignment for a small number of qubits. This instances are 
produced by adding clauses at random until there is exactly 
only one satisfying assignment. In order to apply adiabatic 
quantum computation the initial Hamiltonian Hq taken is a 
magnetic field in the a; direction 



Ho 



(72) 



i=l 



where di is the number of clauses in which qubit i appears. 
Then, for each instance, the ground state is computed for sev- 
eral values of s of the Hamiltonian, H{s) = (1 — s)Ho + sHp 
and its corresponding entanglement entropy of half a chain. 
The mean of the entanglement entropy over these 300 in- 
stances is performed and plotted respect to the s parameter 
for different sizes of the system in Fig. |6| We can observe a 
peak of the entropy around the critical value Sc ^ 0.7. 

We interpret this behaviour of the entanglement entropy as 
follows: initially the system is in a product state and its en- 
tanglement is zero. Then, the evolution makes the system ex- 
plore different solutions by means of superposition states of 
them, that is, it becomes more and more entangled. Finally, 
the system throws away the bad solutions, the entanglement 
decreases, until the best solution is found and it rests in a prod- 



uct state again. Roughly speaking, the power of the quantum 
computer respect to the classical one underlies in the paral- 
lelism during the computation that the superposition principle 
allows. 

Let us make some warning remarks. The numerical simu- 
lations performed for the Exact Cover problem cannot deter- 
mine the complexity class for the quantum algorithm. It is 
generally believed that quantum computers will no be able to 
handle NP-complete problems. Yet, the simulation shows that 
the best this quantum algorithm can achieved still requires a 
huge amount of entanglement in the register. 

The divergence of the entanglement entropy that occurs at 
the critical point Sc ha ve als o been observed in Shor's factor- 
ing algorithm in Ref. ll ISll . where entropy grows exponen- 
tially fast respect the number of qudits. This, again, makes 
this algorithm hard to simulate classically. 

Notice that in the solution of other problems, the explosion 
of the entropy could occur at Sc = 1 in such a way that the en- 
tanglement entropy was monotonically increasing. This sim- 
ilar behaviour of the entropy to the quantum phase transition 
is, therefore, problem dependent. 

Recently, there has been appeared a new quantum algorithm 
for SAT problems that improve the previous results. It consists 
of a hybrid procedu re tha t alternates non-adiabatic evolution 
with adiabatic steps ll 19ll . 



C. One way quantum computation 

The one-way quantum computation (or measurement based 
QC) is a method to perform quantum computation that con- 
sists of: (/) first, an entangled resource state is prepared, and 
(//) then single qubit measurements are performed on it. It is 
called "one-way" because the entanglement of the state, which 
is the resource of the quantum computation, is destroyed by 
the measurements as the computation is being performed. Al- 
though the output of each individual measurement is random, 
they are related in such a way that the computation always 
succeeds. The idea is that depending on the previous outcome, 
one chooses the basis for the next measurements. This implies 
that the measurements cannot be performed at the same time. 

This kind of computation was introduced in Ref. Il20ll 
where there was shown that with an initial particular state, 
called cluster state, any quantum computation could be simu- 
lated. Later on, other useful s tates to perform one-way quan- 
tum computation were found lfl2ll [m [l2l [1241 \l2S\ 

The fact that the measurement based quantum computation 
is universal is non-obvious, since a quantum computation is 
a unitary process, while a measurement is a random process. 
The key point is that there are two kinds of qubits in the spin 
system: the cluster qubits which will be measured in the pro- 
cess of computation and the logical qubits which constitute 
the quantum information that is going to be processed. 

Although, globally, entanglement is expected to decrease 
along the quantum computation due to the single qubit mea- 
surements, in the set of logical qubits (the register that will be 
read out at the end of the computation), the entanglement may 
increase. Notice that if the initial state fulfils an area law, the 
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entanglement is enough that the register of the logical qubits 
is as entangled as possible. That is, area law on a 2D state is 
just what is needed to have linear maximal entanglement on a 
register defined on a line in that state. Cluster states are just 
enough to handle the expected entanglement in the register 
In this respect, it has been recently shown that most quantum 
states are too entangled to be usefu l in o rder to perform mea- 
sured based quantum computation lll26ll . 

Vm. CONCLUSION: ENTANGLEMENT AS THE 
BARRIER FOR CLASSICAL SIMULATIONS 

Entanglement is the genuine quantum property that escapes 
classical physics. The Hilbert space structure of a multi- 
partite quantum system allows for superpositions of exponen- 
tially many elements of the basis. Entropy of entanglement is 
a way to quantify the amount of quantum correlation between 
parts of such a multi-partite system. Entanglement entropy 
is, then, a genuine measure of the global quantumness of the 
state. 

It serves as a conclusion to recall the deep implications of 
entanglement entropy in the possibility of producing faithful 
classical simulations of quantum mechanics. In Ref. IIIO6II . 
it was proven that efficient simulations are possible for any 
system where all its Schmidt decompositions in two arbitrary 
parts would carry little entropy. Therefore, entanglement is at 
the heart of the separation between efficient and non-efficient 
simulations of quantum mechanics. 

What is not fully understood is what is the best strategy 
to classically account for quantum correlations. Two general 



and clever ideas are available in the literature. The first idea 
consists in exploiting the fact that typical interactions are lo- 
cal. This suggests that entanglement should be created se- 
quentially in space from each local degree of freedom to its 
nearest neighbours. Then, a one-dimensional state can be rep- 
resen ted as a matrix product state which captures such a prin- 
ciple I127II . In higher dimensions, sta tes ca n be represented as 
Projected Entangled Pairs (see Refs. Il28ll ). The second idea 
to classically represent quantum states as efficiently as pos- 
sible consists in reconstructing the correlations in the system 
as a renormalization group tree. This goes under the name of 
Multiscale Entanglement Renormalization Ansatz (MERA) is 
a more sophisticated representation which is specially suited 
for critical systems. The accuracy of the approximations can 
be quantify using the amount of entro py of entanglement that 
the approximation can accommodate ll29ll . 

Multi-partite entanglement branches in many others sub- 
jects that escape this short review. Very likely, much more 
work is still needed to get a profound understanding of the 
role of entanglement in highly structured quantum systems. 
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